B34  J.  Opt.  Soc.  Am.  B/Vol.  26,  No.  12/December  2009 


Ku  et  al. 


Bianisotropic  negative-index  metamaterial 
embedded  in  a  symmetric  medium 


Zahyun  Ku,1  Keshav  M.  Dani,2  Prashanth  C.  Upadhya,2  and  S.  R.  J.  Brueck1’* 

1  Center  for  High  Technology  Materials  and  Electrical  and  Computer  Engineering  Department, 
University  of  New  Mexico,  Albuquerque,  New  Mexico  87106,  USA 
2 Center  for  Integrated  Nanotechnologies,  Los  Alamos  National  Laboratory,  Los  Alamos, 

New  Mexico  87545,  USA 
Corresponding  author:  brueck@chtm.unm.edu 


Received  August  3,  2009;  revised  August  5,  2009;  accepted  August  12,  2009; 
posted  August  13,  2009  (Doc.  ID  115231);  published  September  4,  2009 

In  order  to  more  clearly  observe  the  bianisotropic  effects  due  to  fabrication-induced  structural  asymmetries  in 
negative-index  metamaterials  based  on  a  fishnet  structure,  it  is  necessary  to  measure  the  optical  properties 
with  symmetric  substrate  and  superstrate  bounding  layers.  This  is  accomplished  in  this  report  using  an  index¬ 
matching  fluid  and  identical  substrate  and  superstrate  glass  materials.  ©  2009  Optical  Society  of  America 
OCIS  codes:  160.3918,  220.4241,  350.3618,  350.4238. 


1.  INTRODUCTION 

For  many  applications,  thick  (multiple-functional  layer) 
low-loss  double  negative-index  metamaterials  (NIMs)  sat¬ 
isfying  the  sufficient  condition  (Re(e)  =  e'<0;  Re(pi)  =  pi' 
<0),  where  s  (pi)  is  the  permittivity  (permeability),  are 
desirable  [1].  To  achieve  optical-frequency  NIMs,  the  as¬ 
pect  ratio  (height/width)  needs  to  be  large.  For  instance, 
in  previous  reports  the  NIMs  based  on  the  fishnet  struc¬ 
ture  with  negative  refractive  index  at  near-infrared  and 
visible  frequency  were  fabricated  with  aspect  ratios  of 
—0.4  and  —1.4,  respectively  [2,3].  Present  fabrication  of 
fishnet  NIMs  largely  employs  standard  semiconductor 
techniques:  conventional  lithography  (interferometric, 
e-beam,  focused-ion-beam,  or  nano-imprint),  electron- 
beam  evaporation,  etching,  and  liftoff  processing  [2-7]. 
Several  reports  have  shown  that  a  nonnegligible  sidewall 
angle  results  from  this  process  sequence  [8-10] .  This  side- 
wall  angle  results  in  a  bianisotropic  optical  response 
(with  different  reflectivities  from  the  substrate  and  super¬ 
strate  directions).  With  the  exception  of  the  focused-ion- 
beam  processing  [5] ,  all  reported  fishnet  NIMs  have  been 
fabricated  with  a  lift-off  process  and  sequential  deposition 
of  metal  and  dielectric  films.  Because  materials  accumu¬ 
late  on  the  tops  of  the  posts  that  define  the  structure, 
each  successive  film  deposition  has  a  somewhat  larger  ap¬ 
erture  on  the  bottom  metamaterial  film,  giving  rise  to  a 
nonzero  sidewall  angle  and  to  optical  bianisotropy  as 
shown  in  Figs.  1(a)  and  1(b).  This  sidewall  angle  will  vary 
depending  on  the  height  of  the  sacrificial  posts  and  on  the 
deposition  geometry  (distance  from  the  source  and  colli- 
mation).  Etching  approaches  to  forming  the  metamaterial 
are  also  subject  to  top-to-bottom  variations  of  the  film 
stack  geometry. 

The  negative  permeability  in  fishnet  structures  is  as¬ 
cribed  to  the  magnetic  response  from  a  cut-wire  portion  of 


the  structure  and  the  negative  permittivity  to  the  electric 
response  in  the  thin-metal  wire  portion  of  the  fishnet 
NIM,  which  together  lead  to  a  negative  refractive  index 
[11,12].  However,  in  a  structure  with  a  non-negligible 
sidewall  angle,  the  cut-wire  structure  magnetic  response 
induced  by  H  is  coupled  with  an  electric  response,  and, 
similarly,  in  thin-metal  wires  the  electric  response  in¬ 
duced  by  E  is  coupled  with  a  magnetic  response,  which 
gives  rise  to  bianisotropy  and  to  differences  in  the  reflec¬ 
tivity  measured  from  the  substrate  and  superstrate  sides 
of  the  metamaterial  [8,13-18].  We  have  previously  re¬ 
ported  on  this  bianisotropic  behavior  for  an  asymmetric 
case  of  a  glass  substrate  and  an  air  superstrate  [8].  Here 
we  extend  these  measurements  and  simulations  to  the 
symmetric  case  of  a  glass  substrate,  an  index-matching 
fluid  (refractive  index-matching  liquids,  Cargille  Lab, 
Inc.)  in  the  metamaterial  spaces,  and  a  glass  superstrate 
as  shown  in  Fig.  1(c).  This  symmetric  cladding  approach 
provides  a  clearer  observation  of  the  bianisotropic  charac¬ 
teristics  resulting  from  the  nonzero  sidewall  angle. 

2.  DESIGN  AND  FABRICATION 

The  structure  of  all  samples  reported  here  consists  of  BK7 
glass  substrate  with  one,  two,  or  three  metal-dielectric- 
metal  functional  layers  with  an  interpenetrating  two- 
dimensional  square  array  of  elliptical  apertures.  The  geo¬ 
metrical  parameters  of  multifunctional  layered  elliptical 
NIMs  are  indicated  in  Fig.  1(c).  The  orthogonal  pitches  of 
the  two-dimensional  gratings  are  both  fixed  at  810  nm 
( p ).  The  thickness  of  the  gold  and  aluminum  oxide  films 
are  fixed  at  25  and  58  nm,  respectively.  The  linewidths  of 
the  gratings  along  x,  y  axes  are  458  and  274  nm,  respec¬ 
tively.  That  is,  the  elliptical  aperture  size  (2a, 26)  is  fixed 
at  536  and  352  nm  (—0.66  p,  —0.44  p  of  the  pitch  for  the 
major/minor  axis). 
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Fig.  1.  (Color  online)  Schematic  view  of  sidewall-angle  formation  mechanism,  (a)  Illustration  of  tendency  for  tops  of  polymer  posts  to  be 
larger  in  e-beam  evaporation  of  constituent  materials,  (b)  Final  structure  with  sidewall  angle  by  liftoff  processing,  (c)  Dimensions  of 
two-functional-layered  superstrate-elliptical  NIM:  p  =  810  nm,  2a  =  536  nm,  26  =  352  nm,  and  6  (sidewall  angle)  =  13°.  The  polarization  di¬ 
rection  is  described  (E,  H,  and  k  denote  electric,  magnetic  field,  and  propagation  direction,  respectively). 


Detailed  fabrication  steps  for  one-  to  three-functional- 
layered  elliptical  negative-index  metamaterials  were: 

•  2.5  X  2.5  cm2,  1-mm-thick  BK7  glass  was  cleaned 
with  piranha  solution  (4:1  volume  concentration  ratio  of 
H2S04  and  H202)  to  remove  any  residual  contamination 
for  10  min,  followed  by  dehydration  processing  on  150  °C 
hotplate  for  6  min. 

•  A  bottom  antireflection  coating  layer  (BARC, 
—  160-nm-thick,  XHRIC-16,  Brewer  Science,  Inc.)  for 
i-line  lithography  was  spun  onto  a  substrate  at  4K  rpm 
for  30  s  and  oven  baked  at  175  °C  for  3  min;  this  was  re¬ 
peated  three  times  to  make  tall  posts;  negative-tone  pho¬ 
toresist  (NR7-500P,  Futurerex)  was  spread  atop  BARC 
layers  at  4K  rpm  for  30  s  and  hotplate-baked  at  150  °C  for 
60s. 

•  Interferometric  lithography  using  a  355  nm  third- 
harmonic  YAG:Nd3+  laser  source  produced  a  periodic  el¬ 
liptical  hole  pattern  in  the  photoresist  layer  with  two  suc¬ 
cessive  orthogonal  one-dimensional  exposures  with  power 
50  mJ  and  different  doses  for  asymmetrical  aperture 
shape;  a  2  min  post-bake  at  110°C  was  used;  the  exposed 
and  post-baked  sample  was  developed  (MF-702,  Shipley, 
Inc.)  for  13  s  and  rinsed  with  deionized  water  as  shown  in 
Fig.  2(a). 


Fig.  2.  Scanning  electron  microscope  images  of  (a)  a  periodic  el¬ 
liptical  hole  PR  pattern  defined  by  interferometric  lithography, 
(b),  (c)  BARC  posts  etched  by  anisotropic  02  plasma  RIE  using  Cr 
as  the  selective  etching  mask,  (d)  25-nm-thick  Au  deposition  af¬ 
ter  (b),  (c). 


•  After  the  PR  pattern  was  defined,  25  nm  Cr  was  de¬ 
posited  with  an  e-beam  evaporator  at  ~10-7  mTorr  and 
—0.05  nm/s  as  a  selective  etching  mask,  followed  by  lift¬ 
off  processing  with  acetone  to  remove  the  PR  layer.  This 
resulted  in  an  array  of  elliptical  Cr  disks  atop  the  BARC 
layer. 

•  An  anisotropic  02  plasma-reactive  ion  etch  was  used 
to  transfer  the  Cr  pattern  into  the  BARC  layer  with  a 
small  isotropic  etching  component  to  provide  a  slight  un¬ 
dercut  for  the  final  lift-off  step  as  shown  in  Figs.  2(b)  and 
2(c). 

•  Au  (25  nm)/Al203  (58  nm)  were  deposited  alterna¬ 
tively  by  e-beam  evaporation  for  one  to  three  functional 
layers  as  shown  in  Fig.  2(d). 

•  Finally,  an  02  plasma  ash  was  used  to  remove  the 
BARC  posts,  which  leads  to  the  final  structure  with  thick¬ 
ness  of  one  to  three  functional  layers,  108  nm,  191  nm 
and  274  nm,  respectively.  The  insets  of  Fig.  3  show  the  fi¬ 
nal  multifunctional  layered  elliptical  negative-index 
metamaterial  structures. 


Wavelength  (pm) 


Fig.  3.  (Color  online)  Measured/simulated  transmission,  reflec¬ 
tance  (from  superstrate  and  substrate),  and  absorption  (1-T-R) 
for  superstrate-elliptical  negative-index  metamaterials,  depend¬ 
ing  on  the  propagation  direction,  with  the  polarization  shown  in 
Fig.  1(c).  Shaded  areas  represent  wavelength  regions  of  negative 
refractive  index  from  simulation  (see  below).  Insets  display  SEM 
images  with  one-  to  three-functional-layered  eNIMs. 
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3.  FTIR  MEASUREMENT  AND  SIMULATION 

The  transmission  and  reflectance  from  the  superstrate 
and  substrate  sides  were  measured  with  a  Nicolet  Fourier 
transform  infrared  spectrometer  (FTIR)  with  the  specific 
polarization  shown  in  Fig.  1(c),  a  quartz  beam  splitter, 
and  a  DTGS-KBr  detector.  The  transmission  was  mea¬ 
sured  at  normal  incidence  and  normalized  to  the  spec¬ 
trum  of  a  bare  glass  substrate.  The  reflectance  measure¬ 
ment  was  carried  out  with  an  11°  off-normal  incident 
beam  and  was  normalized  to  the  spectrum  reflected  from 
a  gold  mirror.  The  first  and  third  rows  of  Fig.  3  show  mea¬ 
sured  transmission,  reflectance  and,  absorption  (1-T-R) 
depending  on  the  propagation  direction  (superstrate  to 
substrate  and  vice  versa)  of  superstrate-elliptical 
negative-index  metamaterial  (s-eNIM)  with  one  to  three 
functional  layers. 

Commercial  software  based  on  a  finite  integration  tech¬ 
nique  (CST  Microwave  studio:  Computer  Simulation 
Technology  GmbH,  Darmstadt,  Germany)  was  used  to  ob¬ 
tain  the  electromagnetic  response  of  one-  to  three- 
functional-layered  s-eNIMs  depending  on  the  propagation 
direction  [19].  We  used  perfect  electric  conductor  (PEC) 
and  perfect  magnetic  conductor  (PMC)  boundary  condi¬ 
tions  between  unit  cells  with  ttsubstrate  =  7*  superstrate 
~ ^index-match  liquid  ~  1*5,  ftAi203=L62  and  the  Drude  model 
for  gold  permittivity,  where  o)p  (plasma  frequency) 
=  9.02  eV,  coc  (collision  frequency)  =  0.081  eV.  This  colli¬ 
sion  frequency  is  increased  by  a  factor  of  three  compared 
to  that  of  bulk  gold  to  account  for  additional  scattering 
mechanisms  in  this  polycrystalline  thin  film  structure 
[20] .  As  a  result  of  the  broken  structural  symmetry  asso¬ 
ciated  with  the  nonzero  sidewall  angle  (Rsuperstrate 
ARSubstrate)>  modified  retrieval  methods  [8,13-18]  should 
be  used  to  extract  the  effective  parameters  rather  than 
the  standard  method  associated  with  the  symmetric 
structure  along  the  propagation  direction  [21]. 

For  this  defined  polarization,  the  electromagnetic  re¬ 
sponse  of  sidewall-angled  NIMs  can  be  written  as  a  one- 


Fig.  4.  (Color  online)  n'  ,s'  ,jul'  necessary  condition  (s' •  jul" 
+  jul’-s"<  0),  and  sufficient  condition  (s' ,  jul' <  0)  for  negative  re¬ 
fractive  index  (ft'<0)  for  one-  to  three-functional-layered 
s-eNIMs  with  sidewall  angles  (4°,  13°,  22°)  are  plotted.  has 
been  multiplied  by  10  for  clarity. 


dimensional  constitutive  equation  resulting  in  D  =  s0s-E 
-(i/c)£-H,  B  =  (i/c)£;-E+ julqjul-H,  where  c,  e0,  and  julq  are 
the  speed  of  light,  permittivity,  and  permeability  of  free 
space.  The  material  bianisotropy  parameter  £  is  related 
the  constitutive  parameters  by  n2  =  sQs-  julq/jl-1;2  [8,13-18]. 
The  good  agreement  between  experiment  and  simulation 
as  shown  in  Fig.  3  gives  confidence  to  investigating  the  bi¬ 
anisotropy  effects  through  effective  parameters. 

In  Fig.  4,  we  present  the  real  parts  of  the  effective  pa¬ 
rameters  (n' ,s' ,jul'),  the  real  part  of  the  material  bi¬ 
anisotropy  parameter  (£'),  the  necessary  (s' •  jul"  +  jul' •  s" 
<0)  and  sufficient  conditions  (s' ,  jul'  <  0)  for  a  negative  re¬ 
fractive  index  for  one-  to  three-functional-layered 
s-eNIMs  with  sidewall  angles  of  4°;  13°;  and  22°.  By  cau¬ 
sality  (n">0),  the  negative  refractive  index  region  (n' 
<0)  is  equivalent  to  the  region  of  s'  •  jul"  +  jul'  - s" <  0.  The 
bandwidth  of  the  sufficient  condition  (shaded  region)  as¬ 
sociated  with  a  double-negative  material  or  a  negative  in¬ 
dex  with  lower  loss  becomes  narrower  as  the  sidewall- 
angle  increases.  Specifically,  the  region  of  the  sufficient 
condition  is  limited  by  the  bandwidth  of  the  negative  per¬ 
meability  (jh'  <  0)  in  the  case  of  a  sidewall-angled  one- 
functional-layered  s-eNIM  and  a  two-functional-layered 
s-eNIM  with  sidewall  angles  (4°  and  13°),  as  shown  in 
Figs.  4(a)— 4(e)  and  it  is  influenced  by  the  region  of  nega¬ 
tive  permittivity  (e'<0)  in  the  case  of  a  two-functional- 
layered  s-eNIM  with  sidewall  angle  (22°)  and  all  sidewall¬ 
angled  three-functional-layered  s-eNIMs  as  shown  in 
Figs.  4(f)-4(i). 

In  general,  the  double-negative  bandwidth  depends  pri¬ 
marily  on  the  bandwidth  of  the  resonant  negative  perme¬ 
ability;  however,  s-eNIMs  with  both  two  and  three  func¬ 
tional  layers  and  larger  sidewall  angles  as  shown  in  Figs. 
4(f)-4(h)  are  impacted  well  by  the  vertical  shift  of  the 
negative  effective  permittivity  due  to  a  shifted  effective 
plasma  frequency  (o)p).  In  Fig.  4(i),  a  22°  sidewall-angled 
three-functional-layered  s-eNIM  has  no  region  of  double¬ 
negative  material  because  the  regions  of  negative  perme¬ 
ability  and  permittivity  no  longer  overlap.  The  necessary 
condition  is  also  not  satisfied,  so  the  region  of  negative  in¬ 
dex  is  completely  eliminated  by  the  effects  of  the  sidewall 
angle.  In  Fig.  4,  it  is  obvious  that  effective  plasma  fre¬ 
quency  (o)p)  tends  to  decrease  as  the  sidewall  angle  in¬ 
creases  because  it  is  related  to  the  fraction  of  space  occu¬ 
pied  by  the  metal  plates  [22] .  Also,  the  antiresonant  peak 
of  the  effective  permittivity  becomes  weaker  as  the  side- 
wall  angle  increases  because  the  electric  response  be¬ 
comes  weaker,  which  is  a  consequence  of  the  unequal  size 
of  the  metal  plates.  The  antiresonant  behavior  of  the  real 
part  of  £  becomes  stronger  with  an  increase  of  the  side- 
wall  angle.  At  the  magnetic  resonance  for  a  zero  sidewall¬ 
angled  (ideal)  cut-wire  structure,  the  magnetic  field  in¬ 
duces  an  antisymmetric  current  loop  as  shown  in  Figs. 
5(a)  and  6(a),  producing  a  B  field  confined  mainly  between 
the  metal  plates  as  shown  in  Figs.  5(b)  and  6(b).  Figures 
5(c)  and  6(c)  show  that  the  electric  field,  resulting  from 
accumulated  charges  at  the  edges  of  each  metal  plate  in¬ 
duced  by  the  antisymmetric  currents,  is  expected  to  be 
antisymmetric  with  respect  to  the  y  =  0  plane,  to  be  sym¬ 
metric  with  respect  to  the  ry-plane  through  middle  of  the 
structure,  and  to  be  confined  between  and  near  the  ends 
of  the  metal  plates.  In  a  zero  sidewall-angle  thin-metal 
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Fig.  5.  (a)  Magnitude  of  y  component  of  current  density  (|e7y|);  (b)  x  component  of  magnetic  field  (\HX\);  and  (c)  z  component  of  electric 

field  (\EZ\)  at  a  frequency  of  min{Re[/u,eff]}  for  one-functional-layer  s-eNIM  with  sidewall  angles  (0°,  22°)  are  plotted  along  lines  at  the  edge 
of  cut-wire  structure  (x  =  405  nm).  (d)  \Jy\  along  the  lines  in  edge  of  thin-metal  wire  structure  (y  =  405  nm).  The  values  of  \Jy\,  \HX\,  and  \EZ\ 
have  been  divided  by  1012  [A/m2],  105  [A/m],  and  108  [V/m],  respectively. 


wire  structures  at  electric  resonance,  the  induced  cur¬ 
rents  in  the  upper  and  lower  metal  plates  generated  by  E 
are  parallel,  so  the  oppositely  directed  induced  magnetic 
field  lines  go  around  the  metal  plates  and  cancel  within 
the  space  between  the  metal  plates,  as  shown  in  Figs.  5(d) 
and  6(d).  On  the  other  hand,  for  nonzero  sidewall-angle, 
cut-wire,  and  thin-metal  wire  structures,  the  larger  the 
size  disparity  (increase  of  sidewall  angle),  the  stronger 
the  net  electric  dipole  response  in  the  cut-wire  structure 
and  the  magnetic  dipole  response  induced  by  asymmetric 
currents  in  the  thin-metal  wires.  The  reason  is  that  for 
the  nonzero  sidewall-angled  cut-wire  part,  the  magnitude 


of  the  vector  summation  of  electric  dipoles  tends  to  be 
larger  as  the  fabrication-induced  sidewall  angle  in¬ 
creases,  resulting  from  unbalanced  current  density  distri¬ 
bution  induced  by  H.  In  other  words,  there  is  no  longer  an 
antisymmetric  current  density  with  respect  to  the  xy 
plane  through  the  middle  of  the  structure,  as  shown  in 
Figs.  5(a)  and  6(a).  This  also  results  in  an  asymmetric 
electric  field  with  respect  to  the  xy  plane  through  the 
middle  of  the  structure  because  of  the  accumulated 
charge  distribution  at  each  metal  plate,  as  shown  in  Figs. 
5(c)  and  6(c).  For  the  nonzero  sidewall-angled  thin-metal 
wires  part,  the  magnitude  of  the  induced  magnetic  dipole 
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Fig.  6.  (a)  Magnitude  ofy  component  of  current  density  (|Jy|);  (b)  x  component  of  magnetic  field  (\HX\);  and  (c)  z  component  of  electric 

field  (| Ez\)  at  a  frequency  of  minIRef/cteff]}  for  two-functional-layered  s-eNIM  with  sidewall  angles  (0°,  22°)  are  plotted  along  lines  at  the 
edge  of  cut-wire  structure  (x  =  405  nm).  (d)  \Jy\  along  the  lines  in  edge  of  thin-metal  wires  structure  (y  =  405  nm).  The  values  of  \Jy\,  \HX\, 
and  \EZ\  have  been  divided  by  1012  [A/m2],  105  [A/m],  and  108  [V/m],  respectively. 
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moment  between  the  metal  plates  tends  to  be  stronger 
owing  to  a  larger  differential  strength  of  induced  currents 
by  E  in  the  upper  and  lower  metal  plates,  as  shown  in 
Figs.  5(d)  and  6(d).  Notice  that  in  addition  to  the  depen¬ 
dence  on  the  sidewall  angle,  the  antiresonant  behavior  of 
Re(£)  with  the  same  sidewall  angle  as  shown  in  Fig.  4  also 
becomes  stronger  as  the  number  of  functional  layers  in¬ 
creases  [17]. 


4.  SUMMARY 

In  summary,  we  have  experimentally  demonstrated  the 
optical  effects  of  fabrication-induced  structural  asymme¬ 
try  (sidewall  angle)  in  one-  to  three-functional-layered  el¬ 
liptical  negative-index  metamaterials  embedded  in  a  sym¬ 
metric  medium  consisting  of  a  glass  substrate,  a  glass 
superstrate,  and  an  index-matching  liquid.  Numerical 
simulations  are  in  good  agreement  with  the  experiment. 
For  application,  negative-index  metamaterials  should  be 
improved  through  minimization  of  the  structural  asym¬ 
metry  to  achieve  the  double-negative  material  or  multi¬ 
functional  layered  structures  (the  proposed  fabrication 
procedure  was  reported  in  [23]). 
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